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Abstract 

We study the quantum phase transition occurring in an infinite homogeneous sys- 
tem of spin 1/2 fermions in a non-relativistic context. As an example we consider 
neutrons interacting through a simple spin-spin Heisenberg force. The two critical 
values of the coupling strength — signaling the onset into the system of a finite 
magnetization and of the total magnetization, respectively — are found and their 
dependence upon the range of the interaction is explored. The spin response func- 
tion of the system in the region where the spin-rotational symmetry is spontaneously 
broken is also studied. For a ferromagnetic interaction the spin response along the 
direction of the spontaneous magnetization occurs in the particle-hole continuum 
and displays, for not too large momentum transfers, two distinct peaks. The re- 
sponse along the direction orthogonal to the spontaneous magnetization displays 
instead, beyond a softened and depleted particle-hole continuum, a collective mode 
to be identified with a Goldstone boson of type II. Notably, the random phase ap- 
proximation on a Hartree-Fock basis accounts for it, in particular for its quadratic 
— close to the origin — dispersion relation. It is shown that the Goldstone boson 
contributes to the saturation of the energy-weighted sum rule for « 25% when the 
system becomes fully magnetized (that is in correspondence of the upper critical 
value of the interaction strength) and continues to grow as the interaction strength 
increases. 

Key words: spontaneous symmetry breaking, response functions, random phase 
approximation 
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1 Introduction 



As emphasized by lachello [1], atomic nuclei offer an exciting ground to explore 
phase transitions occurring at different temperatures T. Indeed, for these sys- 
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terns at T = a phase transition is signaled by the change they undergo from 
spherical to ellipsoidal shapes; at T ^ 20 MeV evidence exists for a liquid-gas 
transition and, finally, at T pa 200 MeV the deconfinement of the hadronic 
(nuclear) matter is conjectured to take place. 

As it is well-known, at T = the phase transitions are referred to as quantum 
phase transitions. In fact, they occur in correspondence to special (critical) 
values of a coupling constant, viewed as a continuous variable, entering into 
the nuclear Hamiltonian: thus, they are driven, rather than by thermal fluc- 
tuations, as it is the case at finite T, by quantum fluctuations. 

In connection with quantum fluctuations, besides atomic nuclei, also neutron 
stars are systems worth to be explored. Indeed, the strong magnetic field 
present in these systems might be generated by a strong ferromagnetic core in 
the stellar interior. Although a lot of effort has been devoted to this issue (see 
Refs. [2,3,4,5,6] for some recent work), no firm conclusions have been reached 
at present. 

In this paper we do not attempt a realistic calculation of the spin polarization 

of a neutron star, but rather we address general aspects of a quantum phase 
transition in an infinite homogeneous system of interacting neutrons, viewing 
the strength Vi of a spin-spin ferromagnetic neutron-neutron force (admittedly 
schematic with no pretense of being realistic) as a control parameter. In Ref. [7] 
we have indeed found that a critical value of the control parameter, namely 
yiower^ cxists such that, for Vi > K/""^"", an incipient ferromagnetism sets in 
into the system. We have computed Vj°'"'^^, following two different paths: 

• Random Phase Approximation (RPA) on a Hartree-Fock (HF) basis (RPA- 
HF), 

• anomalous single-particle propagator formalism, 

the latter only in the case of a zero-range force, finding that they lead to 
identical results. 

As it is well-known, the onset of the ferromagnetic phase into an infinite, 
homogeneous system is signaled, on the one hand, by the divergence of the 
system spin response function at vanishing frequency in the limit {cu — 0, q ^ 
0), corresponding, in coordinate space, to a very large spin-spin correlation 
length at the critical point. On the other hand, in the broken phase the order 
parameter < M > (M being the system's magnetization operator) acquires a 
finite expectation value: in Ref. [7] we obtained the latter through the Stoner 
equation [8]. Equivalently one says that a vanishing (or a finite) value of 
< M > corresponds to a symmetric (or to a broken) vacuum, respectively. 

A second critical value Vi^^^^ of the control parameter, associated to the 
system's full magnetization or to a completely broken vacuum, exists and was 



2 



found in Ref. [7] in the anomalous propagator framework only for the case 
of a zero-range interaction. In the present work we first explore the impact 
of the finite range of the force, embedded in a parameter A, on ViJ''^^'^ . This 
task was already carried out in Ref. [7] for V^"""^^, which was computed as a 
function of X/kp, kp being the Fermi momentum of the system, in the RPA- 
HF framework. Here we employ the anomalous propagator formalism which 
allows not only to test the previous RPA-HF results for Vj°'"'^'^{X/kF), but also 
to get the range dependence of the upper critical value. 

A second issue addressed in this work relates to the nature of the collective 
modes of the system in the broken vacuum. In this connection one should dis- 
tinguish between collective modes in the (longitudinal) 2;-direction (assumed 
to coincide with the one along which the system's spontaneous magnetization 
occurs) and in the direction orthogonal to z (transverse). In the latter case 
the collective excitations truly represents the Goldstone modes of the system: 
these inevitably occur in presence of a spontaneous breaking of a continuous 
symmetry (in the present case the spin-rotational one). We shall later address 
the issue of the number of these modes: here it suffices to say that in the 
present case just one Goldstone mode exists. 

We compute the collective modes of the system in the RPA-HF framework, 
using again, for the sake of simplicity, a zero-range ferromagnetic force and 
determining as a preliminary the region in the {cu, q) plane where the response 
of the system occurs: this turns out to be different for the modes developing 
in the 2;-direction and for those expanding in the plane orthogonal to z. 

In the former the region results from the interplay of four parabolas, charac- 
terized by the Fermi momenta kp and kp, for spin up and spin down particles, 
respectively, which identify the broken vacuum. The values of kp and kp are 
set by the strength of the interaction Vi. The four parabolas eventually coa- 
lesce into two when the system becomes fully magnetized. In these kinematical 
domains the response function to a longitudinal probe displays a collective be- 
havior, however embedded in the particle-hole continuum and hence damped. 
Specifically, when the system is partially magnetized, for not too large values 
of q — namely where the Pauli principle is active — the response displays 
two maxima, one of which is strongly softened and enhanced with respect to 
the free case, whereas the other is pushed up close to the upper boundary 
of the response, where it almost coincides with the peak already displayed 
by the free response: this situation is reminiscent of the splitting of the giant 
dipole resonance in deformed nuclei [9], although the physical interpretation 
of the upper peak is different in the two cases (see later). For larger q — where 
the Pauli principle is no longer felt — the two maxima merge into one, still 
somewhat softened with respect to the free case, if q, while large, is not too 
large. 
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Remarkably, for a fully magnetized system — that is for a completely broken 
vacuum — no collective mode in the 2;-direction turns out to be possible: an 
occurrence, however, strictly valid only for a zero-range interaction. 



Turning to the transverse collective modes, when the vacuum is broken by 
a ferromagnetic spin-spin interaction a region free of single-particle response 
opens up in the corner of the {cu, q) plane: here is where the Goldstone boson 
lives. Actually the latter, as we shall see, displays — in the combined RPA 
and HF framework and for not too large momenta — a parabolic (rather 
than linear) dispersion relation: hence, according to the general theory of the 
spontaneous symmetry breaking in the non-relativistic regime [10], it turns 
out to be a Goldstone mode of second kind. It is remarkable that the RPA- 
HF theory accounts for its existence through a parabolic dispersion relation, 
which, however, for larger momenta deviates from the parabolic behavior. 
Indeed, it develops, just before entering into the particle-hole continuum, a 
maximum that entails the vanishing of the group velocity of the propagating 
mode. 



Much information on the collective motion in a many-body system can be 
gained already through the moments of the response function, referred to as 
sum rules. In the final part of this work we shall consider the zeroth (5*0) and 
the first {Si) moments, commonly called the Coulomb and the energy- weighted 
sum rules, respectively. Specifically, we shall address the issue of finding out 
how their value is affected by the amount of symmetry breaking occurring in 
the vacuum, i. e. in the system's ground state. 



The present paper is organized as follows: in Section 2 we deal with the prob- 
lem of determining the critical values of the control parameters and their 
dependence upon the range of the interaction. In Section 3 we explore the 
response of the system in the presence of a spontaneous symmetry breaking. 
Specifically, we analyze the system's response in the direction (referred to as 
longitudinal) along which its spontaneous magnetization occurs. In Section 4 
we study the system's response orthogonal to the direction of the spontaneous 
magnetization (referred to as transverse). Here the Goldstone nature of the 
system collective excited state is thoroughly discussed. Finally, in Section 5, 
the same issue is addressed in the context of the moments of the system's 
response functions (the so-called sum rules). In the concluding section we 
summarize our work and shortly discuss the nature of the quantum phase 
transitions addressed in our research. Also, we shortly discuss their signifi- 
cance for the physics of the atomic nuclei and of the neutron stars. 



4 



2 The critical values of the control peirameter 



We have recalled in the previous Section that unlike the symmetric vacuum, 
which is specified by just one Fermi momentum kp, the broken vacuum of our 

system is characterized by two Fermi momenta kp and kp-. these fix the densi- 
ties of the particles with spin up and spin down, respectively. The equilibrium 
condition for our system at zero temperature reads then: 

where the single-particle energies 



OUr^— (2b) 

embody the exact Hartree-Fock (HF) self-energy^ that, while diagonal, is no 
longer proportional to the unit matrix. Indeed for the interaction (in momen- 
tum space) 

V{k) = Vrar ■ ^2^^ (3) 
the anomalous HF self-energy turns out to read 
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for spin up neutrons and 



S — (^kp,kp,k) — Yl^-^.(^kp, kp, k) (4b) 

^ We remind the reader that, although our interaction is of pure exchange character, 
an Hartree term arises in the self-energy owing to the broken vacuum (see Ref. [7] ) . 



5 



for the spin down ones. 

Then, inserting (4a) and (4b), evaluated at the appropriate value of k, into 
Eq. (1) one gets the equation 
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which, in the case of a zero-range force, reduces to the simple expression 

{ktf - {kpf 



(5) 



^1 = 



m {k-^f - {kpf 



(6) 



Writing in the above the Fermi momenta k'p and kp in terms of the magneti- 
zation M of the system according to 



k^ = [37rV(l±M)]V3 



(7) 



one obtains the Stoner equation for the critical value of the coupling Vi. In 
the above 
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Stt^ Qtt^ Gtt^ 
is the density of the system and, of course, < M < 1. 



(8) 



Solving Eq. (5) for M — 1, which implies kp — 2^/^kp, and setting x — X/kp 
one gets 
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This formula, though cumbersome, becomes transparent for a zero range in- 
teraction {x — > 00), where it reads 



7r2 



._on mkp 2y^ ^ ' 



6 



B 




0.8 




Fig. 1. In panel (A) the behavior of the critical couplings Vi"^"^^^ (dot-dashed) and 
pdower (^gQ^j^^ computed in the anomalous propagator formalism are shown versus 
X = X/kp; also displayed is in RPA-HF (dashed). In panel (B) one sees Vl°^^'' 

in the effective mass approximation both in the anomalous propagator framework 
(solid) and in RPA-HF (dashed). All curves are in units of —Tr^/mkF- 

a value coinciding with the one found in Ref. [7], while for an infinite range 
force {x 0), where it yields 



lim Vi 



upper 



TT 



1 3 



The behavior of the absolute value of Vi^^*^^ is shown in Fig. 1 (panel A). 
Note that the actual strength of the interaction (3) is T^A^: it would display a 
monotonically increasing behavior with A, showing that the larger the range of 
the force, the smaller the resistance of the medium to become fully magnetized. 
It is worth reminding that this behavior is shared by the 0(3) model [11], 
which enjoys the same symmetry of ours and represents a generalization of 
the Ising model [12] . Also remarkable is that for ranges of the force smaller than 
about (2kp)~^ the critical value of the control parameter becomes essentially 
insensitive to A. 

In concluding this Section we quote the lower critical value of the coupling Vi 
as obtained through the Stoner Equation. For this purpose it merely suffices 
to solve Eq. (5) setting < M >— 0. One gets 
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which for a zero-range interaction reduces to 
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a value coinciding with what found in [7] and for an infinite range force yields 
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\imVl°r''^ ^- (14) 

x-^o rnkp 

Formula (12) is displayed in Fig. 1 (panel A) as a function of x. In the same 
figure it is also shown the result obtained in Ref. [7] through the solution 
of the RPA equation in the long wave length, small frequency domain [13]. 
Notably the two curves differ at intermediate value of x. This outcome relates 
to the approximation made in Ref. [7] , namely of using the effective mass as 
an approximation to the self-consistent HF solution. Indeed, if we do the same 
in the present anomalous propagator framework, we find for the neutron self 
energy 

sf+(4, k^, k) = + (15) 

with 
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(17) 



It is gratifying to see in Fig. 1 (panel B) that now, namely in the effective 
mass approximation, the anomalous propagator framework and the RPA-HF 
lead to identical results. 



3 The system's longitudinal response 



Let us assume the system to undergo a spontaneous symmetry breaking ac- 
quiring a magnetization along the z-axis. We wish to explore the system's 
response to a spin-dependent, but not spin-flipping, external probe acting in 
the z direction, namely described by an operator 

= Y.< r(i\(J^\sa > dlpCisa = ^ aj^a,„(-l)^/^"". (18) 

rp sa 

sa 

For sake of simplicity we confine ourselves to assume a ferromagnetic {Vi < 0), 
spin-dependent, zero-range interaction among neutrons, namely 

V{r)^V,a,-a25{r), (19) 

clearly constant in momentum space (see Eq.(3) in A — > oo limit). 
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We shall compute the response to the probe (18) in the RPA-HF framework 
both in a normal and in a broken vacuum. To this end it is first necessary to 
set up the longitudinal HF anomalous polarization propagator in the broken 
vacuum H^J'^. This is easily achieved starting from the anomalous single- 
particle propagator (the label "b" stands for "broken vacuum") 



^^''"(^,^0) 



G'l^j\k,ko)) 



(20) 



where 



is the propagator for a spin up neutron and 



G^^'^fc, ko) = , ~ + , ^^^^ ~ (21b) 
/co — a;- + IT] ko — cOj^ — irj 

is the propagator for a spin down one and 

ut^ — - (22) 

are the single-particle energies (2) for a zero-range force. In the following we 
shall specify the spin indices for G only, dropping all the other ones, and we 
shall always assume the HF approximation for the single-particle propagator. 

One gets then for the HF polarization propagator in the 2;-direction and in 
the broken vacuum, setting K = (/cq, k) and Q = {u!,q), 
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Fig. 2. The response region in a broken vacuum. In the figure kp = 400 MeV/c 
and kp = 237 MeV/c. The heavy (hght) fine marks the domain where the spin 
up (down) neutrons respond to an external probe. The dashed fines represent the 
boundaries of the Paufi blocked regions. 

r dk 

n«^(g) = n° _(g) = J j^,o{\k + q]- kp)e{kp - k) 



X 



1 



1 



k+q 



.(25) 



Note that the HF expressions for n_|__|_ and 11 are identical to the free ones 

in the case of a zero-range interaction. Moreover both their real and imaginary 
part can easily be computed analytically: one clearly obtains the familiar re- 
sults for a symmetric vacuum [14] with kp replaced by kp and kp, respectively. 

Eqs. (23), (24) and (25) allow us to display in Fig. 2 the response region of the 
system, in the frequency u - momentum q plane, when the vacuum is broken. 
We see in the figure that the global response region is actually made up by 
two response domains: one associated with kp, where the particles with spin 
up respond to the external probe and the other associated with kp, where the 
particles with spin down respond to the external probe. 

Turning to the propagator Uf^^'^^'^ in RPA-HF, it obeys the equation 

n = + n™ vn, (26) 

graphically displayed in Fig. 3, where: 
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Fig. 3. The Dyson equation for the RPA anomalous polarization propagator. 
The solution of Eq. (26) turns out to read (see Appendix A) 



j-|-RPA-HF,6 



1 - Kz(n 
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(28) 



od) 



In the above and Vod correspond to the diagonal and off-diagonal particle- 
hole matrix elements of the interaction (19) in spin space: they are given in 
Appendix A. Noteworthy is that Eq. (28) is formally similar to what one 
gets in systems containing two species of particles, e.g. nucleons and A's (see 
Ref. [15]). 

Equation (28) entails a striking consequence, namely that for a fully broken 
vacuum (a fully magnetized system, for example in the positive z-direction) 
no RPA collective mode exists for a zero-range force. Indeed, in this case, 
since fc]^ = then 11° _ = and also, as shown in Appendix A, Ki=0. The 
situation is clearly illustrated in Figs. 4, 5 and 6, where the system's response 
along the z-axis is shown at g = 5, 50 and 500 MeV/c, respectively. Further- 
more, in each figure, the evolution of the system's response with the amount 
of breaking of the vacuum is also displayed. Accordingly, in each figure the 
responses associated to four pairs of values of and /c]^ arc shown: since the 
density of the system is fixed, these are related by Eq. (8). We further observe 
that each choice of {kp,kp) corresponds to a value of the strength of the in- 
teraction Vi given by Eq. (6). In panel A of all figures the enhancement and 
softening of the response in the symmetric vacuum [kp — 338.13 MeV/c), due 
to the attractive ferromagnetic interaction, is clearly apparent. As one moves 
towards an increasingly broken vacuum and for not too large momenta one 
sees the appearance of a second peak in the response at high energy until, for 
a totally broken vacuum, the collectivity completely disappears in accord with 
the argument given above and the free response is recovered. 



In order to understand the frequency behavior of the response at q 
q — 50 MeV/c it helps to keep in mind that 
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Fig. 4. The response of an infinite neutron's system to a 2;-aligned probe for q = 5 
MeV/c. Panel (A) refers to the case of a symmetric vacuum [kp = kp = kp), 
panels (B) and (C) to a partially aligned broken vacuum, panel (D) to a totally 
broken, fully aligned vacuum. Dotted line: HF (free) response, dashed line: ring 
approximation, solid line: RPA-HF. 
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Fig. 5. The same as in Fig. 4 but for g = 50 MeV/c. 



13 



A 



B 



0.01 



0.008 - 



^ 0.006 
a. 

& 0.004 

a 

^ 0.002 







k:=338.13MeV 

F 

kp=338.13MeV 




k =350 MeV 

F 

k =325.36 MeV 

F 




100 200 300 100 200 300 
CO (MeV) (0 (MeV) 

C D 



0.01 



- 0.008 



0.006 



- 0.004 



0.002 







0.01 



0.008 



k;=400 MeV 
k =237.04 MeV 

F 




k>426.01 MeV 
k =15.72 MeV 

F 



0.01 



0.008 



- 0.006 



0.004 



0.002 



100 200 300 100 200 300 
0) (MeV) CO (MeV) 



Fig. 6. As in Fig. 4 but for q = 500 MeV/c. 
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a) the HF (free) response in the broken, but not fully so, vacuum already 
displays two maxima, when the Pauli principle is active (namely for not too 
large q); 

b) the RPA-HF framework conserves the energy-weighted sum rule even when 
the vacuum is broken (this non trivial result will be further discussed later) . 



Actually, the RPA-HF expression of the response reads 
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V = Ren^^Imn"!: + Imn^FRen"!:. (33) 

In the high energy domain where only the spin up neutrons respond to the 
external probe (namely, where Imll^^ = 0) the above become 
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and 
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(35) 



respectively, and since in this regime Rell^^ is a rapidly decreasing function 
of the frequency (see Fig. 7), it follows that the RPA-HF response at large u 
approaches the free one and thus the second peak displayed by the latter still 
shows up. 

By contrast, in ring approximation one has the simpler expressions 



Airing = -—Im(H^^^^ + H!!i) 



and 



Aing = 1 - ^iRe(H«^ + H"!) + V,' [imiUll + H«^) 
which, in the high-frequency regime, become 
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Fig. 7. Reng^ (solid) and -IrnHg^ (dot) as a function of a; at g = 5 McV/c. 
Note that Rell^^ has a maximum at the upper hmit of the corresponding response 
domain. 

and 

i^ring = 1 - V,Re(Ull + n^F) + V.'ilmUll)''. (39) 

Prom the last equation it appears that at large cu the ring response is strongly 
damped. 



4 The system's transverse response 

In this Section we explore the system's response to a probe aligned in the 
direction orthogonal to the axis along which the spontaneous magnetization 
of the system occurs. For definitiveness we choose the probe to act in the 
x-direction. According to the general theory in a non-relativistic context [10], 
we expect here Goldstone modes to show up. Their number should not be less 
than the number of the broken generators of the continuous symmetry, pro- 
vided that the Goldstone bosons of type II are counted twice. In the case we 
are presently investigating, the number of the broken generators is provided 
by the dimensions of the coset 0(3)/0(2), where 0(3) is the rotation group in 
three dimensions. This group leaves invariant the Hamiltonian of our system 
of interacting neutrons, whereas 0(2) is the rotation group in two dimensions 
and represents the surviving symmetry after the spontaneous breaking has 
occurred. Hence in our case two generators are broken. Accordingly this situ- 
ation is compatible with the existence either of two Goldstone bosons of type 
I — characterized by a dispersion relation linear in the momentum — or with 
the existence of one type II Goldstone boson — which has a dispersion relation 
quadratic in the momentum. 
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Fig. 8. The diagrams corresponding to U^+{Q) (A) and U^^(Q) (B). 



As we shall see, the latter is actually the occurring case for our system. This 
is hardly surprising since indeed type II Goldstone bosons are specific of a 
non-relativistic theory as exemplified by the rotational bands of atomic nu- 
clei, whose energy depends quadratically on the angular momentum (nuclei 
are finite systems!) [9] and by the pairing additional and removal modes in su- 
perconducting nuclei [16] , whose energy depends quadratically on the number 
of pairs. Clearly, in the two above mentioned examples the broken symmetries 
are the rotational and the global gauge ones, which are spontaneously broken 
in deformed and superconducting nuclei, respectively. It would be interesting 
to study a suitable relativistic generalization of our system, since in that case 
the previously mentioned counting rule for the number of Goldstone bosons 
no longer holds [17]. 

In searching for this Goldstone bosons we first ask: where do they live? To 
answer this question we need to consider the transverse HF polarization prop- 
agator 



= / ^,[G^4K)G^-{K + Q) + G_4K)G^4K + Q)] 

= n«^(Q) + n^!(g), (40) 

where we have found it convenient to introduce the quantities 11^^ and H^^, 
whose vertices embody the spin operators 

0-± ± i(^y), (41) 

as shown in Fig. 8. In the HF approximation the expressions for H ^ and 

H_|__ are easily deduced starting from the single particle propagators, already 
employed in deducing the response to a longitudinal external probe, given in 
Eqs. (21a) and (21b). One gets: 
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e{\k + q\ - kj-)9{k+ - k) e{kj^ - \k + q\)0{k - k 



UJ + UJ 



+ 



(42a) 



and 



dk 

(2^ 



e{\k + q\ - k+)e{kp - k) e{k+ - \k + q\)e{k - kp] 



(42b) 



From the above formulas, the response region of the infinite, homogeneous 
neutron's system in the (u, q) plane to a spin-flipping probe ((T±) is deduced 
by searching for the region where, e. g., Il^{Q) develops an imaginary part. 
For this purpose we write 



imn^^;(g) = imm+(g) + imni+(Q), 



being 



Imn^+(Q) 



and 



dk 



(43) 



^—e{\k + g1 - k-p)e{k% - k){-n)5{u; + u;i - 



(44a) 



j^,0{k^ -\k + q])9{k - kt){-n)S{u; + ut - uZ^^). 

(44b) 

Hence, the first contribution to the imaginary part of (namely Imll" ^(Q)) 
lives in the domain 



k^q 



k%q 



2m m 2m m 

while the second one (namely Imn^((5)) lives in the domain 

— h Ao; < a; < — \ h Ao;. 

2m m 2m m 



In the above we have set 



Au; 



■27r2 



(45a) 



(45b) 



(46) 
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Fig. 9. The response region associated to 11^^ (heavy hnes) and to (hght 
hnes) . The soUd and dashed lines correspond to the first and second term in which 
we have spht the imaginary part of 11^^, respectively (see Eq. (43)). The values for 
kp and kp are the same as in Fig. 2. 

and in our derivation we have exploited Eq. (6) and the relation 



k 



the last equality being valid for V'^^^'^ < Vi < Vx'^^^'^ . For definitiveness, in 
the following we shall always choose kp > A;^, which implies Aa; > 0. 

The response region of U^^{Q) can be derived along the same lines, yielding, 
instead of Eq. (45), the following expressions: 



il_^_Aa;<u;<^ + ^-Aa; (48a) 
2m m 2m m 

_^_^_Aa;<a;<-^ + ^-Ac.. (48b) 
2m m 2m m 

It is of importance to observe that the response region related to 11^^ (n+^) is 
shifted with respect to the symmetric case upward (downward) by an amount 
Auj that directly reflects the size of the spontaneous breaking of the vacuum. 
The response regions for 11^^^ and n^|^^ are displayed in Fig. 9. 

Concerning the response function, it is remarkable that the following symme- 
try relation holds valid: n^^(g, uj) = 11^^ (g, —uj), as one can see by comparing 
Eqs. (42a) and (42b). Furthermore, note that, at variance with the symmetric 
vacuum case, now, for o; > 0, also the second piece on the right hand side of 
Eqs. (42a) and (42b) contributes to the system's response, the more so, the 
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Fig. 10. The first order direct and exchange terms for n5:+^(Q). 



smaller q is. To complete our analysis we quote in Appendix B the real and 

IF 

-+ 



the imaginary part of 11^, and . 



It should finally be stated that similar results have been obtained in the con- 
text of asymmetric nuclear matter [18,19]. 

We turn now to discuss the RPA equations for n_+ and n+_. The basic 
ingredients required for their deduction are the first order direct and exchange 
diagrams displayed (for n_|_) in Fig. 10. In our case of a zero-range interaction 
one finds 

ngr(g) = 2V,Il%{Q)Illl{Q) (49) 
for the direct term and 

nS^^(Q) = VJlll{Q)\lll{Q) (50) 

for the exchange diagram. Hence the RPA series (which accounts for both 
contributions) can be easily resummed, leading to 



ttRpa-hf/^n _ ngl(Q) 

^^)- i-wji%i{Qy ^^^^ 



To find the dispersion relation of the Goldstone bosons we search for the poles 
(if any) of the expression 

^ [l-3V^inHF(5-',^)][l-3V^inHF(g,a;)] ^^^^ 

for positive real uo. 

From the numerical analysis we have found that of the two factors appearing 
in the denominator of Eq. (52) only the first one (since A;^ > k~^) vanishes for 
just one real and positive value of a; at a given q. In Fig. 11 we display the 
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q (MeV/c) q (MeV/c) q (MeV/c) 

Fig. 11. The dispersion relation of the Goldstone boson for kp = 426.01 MeV/c 
(heavy soHd lines), which corresponds to Vi = —189.25 MeV fm^. Also displayed 
are the response regions: the light solid and dashed lines correspond to Eqs. (45a) 
and (45b), respectively. In panel A one can appreciate how tiny the energy of the 
Goldstone boson is; in panel B, which enlarges panel A, one can assess the domain 
of validity of the parabolic dispersion relation of the Goldstone mode (dot-dashed 
line); in panel C the Goldstone mode is displayed for three different values of the 
interaction strength, namely Vi = -189.25 (solid), -200 (dot) and -300 MeV fm^ 
(dot-dot-dash) . 

solution of the equation 

[1 - SViReU^liq, u)] = 0, (53) 

which we expect to yield the dispersion relation of the Goldstone boson, should 
the RPA be a trustworthy theory for our many-body system. This turns out 
indeed to be the case, since for small q and u the solution of Eq. (53) can be 
analytically expressed through the expansion of Ren^^(g, a;), which reads 



limRen"^(g,a;) = ";^f \ 



+ 



15 (i/ - Ai/)2 



(54) 



where 



and Au 



-mVi 



-^3l 



27r2(A;+)2 



mr - (k-p) 



Prom the above one gets the following dispersion relation, valid for small values 
of q: 



2m* 



(55) 
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with 

-(i-o(i+3e+e)^- ^^^^ 

From Fig. 11 (panels B and C) it appears that the expression (55) actu- 
ally remains valid over a substantial range of momenta. Thus, the solution 
of Eq. (53) truly corresponds to a type II Goldstone boson as it should. Mi- 
croscopically this mode behaves hke a particle-hole excitation (the particle 
being quite heavy). Physically it can be viewed as a twisting of the local spin 
orientation as the collective wave passes through the system [20]. 

Furthermore, and remarkably, it turns out that for Vi > Vi^^'^'' the Goldstone 
mode continues to exist with a dispersion relation that is parabolic over a 
range of momenta becoming larger as Vi increases. For 1^/°"*^^ < Vi < Vi^^'^^ 
the Goldstone mode displays instead an anomalous behavior: in fact, in this 
range of couplings, in correspondence to a specific momentum, the collective 
mode is characterized by a vanishing group velocity. 

It is worth comparing the dispersion relation (55) with the formula for the 
energy levels of a rotational band. They are identical providing one replaces the 
momentum with the quantized angular momentum (of course, both quantities 
should not be too large) and the mass with the moment of inertia. From Fig. 11 
(panel A) it is clearly apparent how tiny the energy of the Goldstone boson is, 
just as it happens for the rotational bands in nuclear and molecular physics. 

In concluding this Section we display in Figs. 12, 13 and 14 the continuum 

response of the system to an x-aligned probe for the same momenta and vacua 
of Figs. 4, 5 and 6, where the response to a 2;-aligned probe was considered. 
We notice that 

i) for a symmetric vacuum the energy-weighted sum rule is patently obeyed; 

ii) the more the vacuum is broken, the more depleted the particle-hole contin- 
uum is; 

iii) in accord with ii) the more the vacuum is broken, the stronger the Gold- 
stone boson becomes. This item will be quantitatively addressed in the next 
section in the sum rule framework. 



5 The moments of the response function 



In this Section we investigate the non-energy-weighted (Sq) and the energy- 
weighted (5"!) sum rules, exploring their behavior when spontaneous symmetry 
breaking occurs in the vacuum. 

Concerning Si, it is well-known [21] that it is given by the following expression 
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Fig. 12. The response of an infinite neutron's system to a x-aligned probe ioi q = 5 
McV/c. Panel (A) refers to the case of a symmetric vacuum (kp = kp = kp), panels 
(B) and (C) to a partially aligned broken vacuum, panel (D) to a totally broken, 
fully aligned vacuum. Dotted line: HF response, solid line: RPA-HF. 
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Fig. 13. The same as in Fig. 12 but for g = 50 MeV/c. 
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Fig. 14. As in Fig. 12 but for q = 500 MeV/c. 
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S^{q)= dujcv ^ = -(0|[O,[i/,O]]|0). (57) 

Jo np 2 

For the density response of a non-interacting gas of fermions of mass m the 

above is indeed fulfilled and yields 



S--{q)^ Tduu^^^^^^^^^f. (58) 
^ Jo np 2m ^ ^ 



In general, however, Eq. (57) is violated by most of the many-body frameworks, 
the remarkable exception being the RPA-HF theory. In fact, the Thouless 
theorem [22] states that if the system's response is computed in RPA-HF and 
the expectation value on the right hand side of Eq. (57) is taken in the HE 
ground state, then the sum rule is fulfilled. We have indeed verified it by 
computing numerically with very good accuracy the left hand side of Eq. (57) 
and by working out the HE expectation value of the double commutator in the 
same equation — which of course yields g^/2m, if one employs the interaction 
(19) and the vertex (18). 

Remarkably, even when the vacuum is broken 5*1 keeps the above value, as it 
can be inferred from the results reported in Tables 1, 2 and 3. Note that this 
outcome could also be proved along the lines followed in Ref. [23] in the case 
of a symmetric vacuum. Conceptually, it relates to the very meaning of Si, 
namely of expressing the particle number conservation (or the global gauge 
invariance) of the theoretical framework and in the non-symmetric vacuum it 
is the rotational — and not the gauge — invariance to be broken. 

In this instance, at variance with the situation where the probe acts in the 
direction of the spontaneous magnetization, when the spin-fiipping probe is di- 
rected orthogonally to the latter, 5*1 is contributed to not only by the particle- 
hole continuum, but by the collective Goldstone mode as well. Indeed, from 
Eq. (51) one has 

jj^nRPA-HF/A^N ^ imng^(Q) 

which, in the region where Imn^^((5) vanishes, yields 



imn 



RPA-HF^ 



TT 




1 
























ui=uj+{q) 



-5[uj - uj+{q)], 



(60) 



with u}+{q) the solution of Eq. (53), that is the Goldstone boson dispersion 
relation. 
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I + 

kp 




Q rt Jri^.- rir 


Sf^ 




Ao; 


338.130028 


0.0111 


(0.0819+0) 


0.013312 


(0.013312+0) 





350 


0.1091 


I'O 00152+0 10754") 


0.9795 


I'O 013295+0 000017'! 


8.860 


400 


0.65549 


(0.000228+0.65527) 


36.248 


(0.012594+0.000718) 


55.279 


426.01 


1.000 


(0.000111+0.99979) 


96.510 


(0.010664+0.002686) 


96.505 


Table 1 












The non-energy-weig 


;hted and energy-weij 


^hted sum 


rules at g = 5 MeV/ 


c, 



corresponding to g^/2m = 0.0133120 MeV. The Fermi momentum = 
338.130028 MeV/c corresponds to a symmetric vacuum, whereas fc^ = 
426.01 MeV/c corresponds to an almost completely broken vacuum (ground state 
fully aligned in spin space). In the columns associated with the RPA-HF theory, 
the first figure represents the contribution to Sq and Si, respectively, arising from 
the particle-hole continuum; the second figure the one arising from the collective 
Goldstone mode. 





cHF 


^RPA-HF 




^RPA-HF 


Aw 


338.130028 


0.1107 


(0.5842+0) 


1.33120 


(1.33120+0) 





350 


0.1374 


(0.5340+0) 


2.29744 


(1.33120+0) 


8.860 


400 


0.6555 


(0.0233+0.6322) 


37.566 


(1.26414+0.06706) 


55.279 


426.01 


0.9999 


(0.0111+0.9888) 


97.827 


(1.07098+0.26025) 


96.505 



Table 2 

The same as in Table 1 but for q = 50 MeV/c, corresponding to q^/2m = 1.33120 
MeV. 



k+ S^^ 5o^PA-HF Sf^ SfPA-HF Aa; 

338.130028 0.907 1.528 133.120 133.120 

350 0.908 1.524 134.087 133.120 8.860 

400 0.951 1.359 169.356 133.120 55.279 

426.01 1.000 1.000 229.616 133.120 96.505 

Table 3 

The same as in Table 1 but for q = 500 MeV/c, corresponding to g^/2m = 133.120 
MeV. 
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Fig. 15. Contribution of the Goldstone mode to the energy- weighted sum rule as a 
function of the interaction strength for g = 5 MeV/c (sohd hne); the curve starts 
at V^°^^^ = —159.16 MeV fm^. The dashed hne represents the saturation value 
q'^/2m = 0.0133120 MeV. 

Actually, this contribution grows with the amount of symmetry breaking in 
the vacuum, which, in turn, grows with the strength of the force Vi. For 
Vi = Vi^'^'^\ namely when the vacuum is fully aligned in spin-space, the 
Goldstone mode accounts for roughly 25% of the energy- weighted sum rule, 
as it can be deduced from the figures reported in Tables 1 and 2. For still larger 
values of Vi, this amounts keeps increasing until, for Vi — > oo, it exhausts the 
sum rule, as it is seen in Fig. 15, where the Goldstone mode contribution to 
5*1 is displayed versus Vi . Of course, for large values of q the Goldstone boson 
no longer exists (Table 3). 

Concerning the non-energy- weighted sum rule in a symmetric vacuum and for 
a non-interacting system of fermions one has the well-known result 



So{q) 



1 - 



12 \kF 



q < 2kF 
q > 2kF, 



(61) 



which is conserved in the HF theory, but not in RPA or RPA-HF. When the 
vacuum is spontaneously broken by our ferromagnetic force, one finds that the 
impact of the Pauli correlations on 5*0 is lowered with respect to the symmetric 
vacuum case and decreases as the amount of the symmetry breaking grows, as 
it is apparent from Tables 1, 2 and 3. In particular, for a fully broken vacuum, 
5*0 is just 1 for any q — that is the value occurring in the symmetric vacuum 
for q > 2kF — both when our system is explored in the longitudinal or in 
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the transverse direction by a spin-dependent probe: the system's constituents 
no longer feel the Pauh principle, as it should be expected since in the fully 
magnetized case only one species of particles is present. 

The reduced influence of the Pauli principle, when the system in only partially 

aligned, with respect to the situation occurring in the symmetric vacuum, can 
be exploited to investigate (using a spin-fiipping probe) how the collectivity 
of the Goldstone mode is affected by the degree of spontaneous symmetry 
breaking of the vacuum. In fact, in Tables 1 and 2 one sees that the less 
effective the Pauh correlations are, the more collective the Goldstone boson 
is. 



6 Conclusions 

In the present study we have dealt with an infinite, non-relativistic, homoge- 
neous system of neutrons interacting through a simple ferromagnetic force of 
Heisenberg type, our aim being to discuss general aspects of the spontaneous 
symmetry breaking associated with a quantum phase transition. Although 
this theme has been much addressed in the past for generic spin 1/2 fermions 
(see, e. g., Ref. [8]), still we felt it useful to analyze in more detail the depen- 
dence upon the interaction range of the critical values of the coupling and the 
excitation spectrum of a system undergoing a quantum phase transition. 

Concerning the latter, we like to remind that the symmetry breaking taking 
place in our system, namely the onset of a permanent, spontaneous magneti- 
zation, stems from a well identified physical source, a situation very different, 
for example, from the electroweak symmetry breaking in the standard model, 
where the responsible for such an occurrence, namely the Higgs boson, is still 
searched for three decades after having been conjectured. Indeed, in our system 
— just like in a ferromagnet the electromagnetic interaction between neigh- 
boring atoms lines up their spins parallel to each other — the force between 
neighboring neutrons produces the same effect, provided that the strength 
of the interaction exceeds some critical value, Vj°^^^. How close the neutrons 
should be for the phase transition to occur or, equivalently, how the range 
of the interaction affects the value of Vi"^'^^7 We have answered this question 
both numerically and analytically, finding out, as expected, that the longer 
the range, the weaker the strength of the coupling should be in order to induce 
the spontaneous breaking of a symmetry. A further natural question to ask is 
how large Vi should be for the system to become fully magnetized. Also for 
this problem we provide a numerical and an analytic answer, again confirming 
the above referred to correlation between strength and range. 

We have next analyzed the spectrum of the system (or, equivalently, the re- 
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sponse functions to external probes acting on the spin of the constituents), 
following its evolution with the amount of spontaneous symmetry breaking 
occurring in the ground state. Probing the system along the direction of the 
spontaneous magnetization, the chief feature of the response function relates 
to the occurrence of two distinct peaks, which reflect both the action of the 
ferromagnetic force and of the Pauli principle. Indeed, the two peaks merge 
at large transferred momenta, where the latter disappears. 

In the direction orthogonal to the magnetization, on the other hand, a collec- 
tive Goldstone mode shows up as required by the theory. Its dispersion rela- 
tion is indeed parabolic, as it should, over a momentum range that increases 
as the coupling strength Vi increases. Actually, for Vi ^ oo, the Goldstone 
boson exhausts the energy-weighted sum rule and its dispersion relation be- 
comes a perfect parabola. Notably, for a strength Vi intermediate between the 
two critical values V}°^^^ and Vi^"^^^, the Goldstone boson dispersion relation 
displays an anomalous non-parabolic behavior, entailing the existence of a 
wavelength associated with a vanishing group velocity of the collective mode. 
We conjecture this to be a distinctive signature of the Goldstone boson in the 
non-relativistic regime. 

In spite of the simplicity of our interaction, which has of course no pretense of 
being realistic, it appears that our research bears significance for the physics 
of the neutron stars, since it explores the extension of the many-body response 
theory to the situation associated with a broken vacuum in a spin space, which 
is required for the assessment of the magnetic field that neutron stars host. 
A lot of work has actually been lately done on this issue: interestingly, it 
appears that simple effective interactions — such as the Skyrme ones — give 
indeed rise to a phase transition of second kind [2,3] , whereas more microscopic 
many-body approaches — such as the Brueckncr-Hartree-Fock formalism [4] or 
quantum simulations [5,6] — give no indication of a quantum phase transition. 
Generally speaking, this striking difference can be related to the different 
predictions these models give for the particle-hole spin interaction at neutron 
star densities: attractive in the Skyrme models and repulsive in calculations 
based on realistic nucleon-nucleon potentials [24]. Unfortunately, at present 
there are no direct phenomenological constraints on this component of the 
effective nuclear interaction at densities relevant for the neutron stars. 

In connection with the possible application of the concept of spontaneous sym- 
metry breaking to the physics of the atomic nuclei, it should be of course real- 
ized that in finite systems the notion of phase transition applies only approx- 
imately. Yet, one is naturally lead to think of heavy nuclei, where the isospin 
symmetry is broken, leading to ground state configurations quite similar to 
the ones we have been considering in the spin space. However, in isospace the 
symmetry is broken explicitly (by the Coulomb force) and not spontaneously 
(as in our case). Thus, although for heavy nuclei the response region of the 
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system in the (u;, g) plane turns out to be very similar in both cases [18], in 
nuclear matter no coUective Goldstone boson should show up in the small uj, 
small q corner of the [u, q) plane. Furthermore, in the isospin particle-hole 
channel — the one appropriate to the present considerations — no attraction 
seems to exist. 



Appendices 



A Probing the system in the symmetric direction: the RPA matrix 
elements 



Here we derive the exact expression for Y\^f^'^. In Eq. (26) we made an ansatz 
on the form of the equation obeyed by the latter. We have now to fix the 
matrix elements of the potential entering into Eq. (27). This can be easily 
done expanding Eq. (26) to first order in Vi and then matching the result 
obtained in this way with the explicit expression coming from the evaluation 
of the first order diagrams in Figs. (A.l) and (A.2). From Eq. (26) it follows 
immediately that 

n = (i-n°t>)-^n'^. (A.i) 

Hence the evaluation of H simply amounts to invert a 2 x 2 matrix. We get: 



(A.2) 



where 



D^i- Vainly + n° _) + n°^n° _(Kf - K'J- (a.3) 

The complete RPA result is thus given by 
j-[RPA ^ ^ + 2H+_ 

= 1 [n% + H" _ + 2n^+no _ ( v;, - y,)] , (a.4) 

which can be expanded to first order in V obtaining 

nw^p^ = n° + + n° _ + n° + n° + 2n^+Kdn° _. (a.s) 

The direct contribution to H^^^^^"^ comes from the diagram in Fig. A.l and 
reads 

n(i)dir = h^+i^iHO + + n° _yin° _ + 2n°+Kdn'L_. (a.6) 
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Fig. A.l. The first order direct term contributing to H^^^. 
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Fig. A. 2. The first order exchange terms contributing to H^^"^. 



The first order exchange contribution comes from the three diagrams in Fig. A. 2 
and reads 



++ ' 



(A.7) 



If one chooses to keep only the direct contribution to the first order result 
for n {ring approximation) then, from the matching with Eq. (A. 5), it follows 
that 



vf' = Vi, v;7 ^Vi ^ v"'' = 



^dir 



dir 



(A.8) 



On the other hand, keeping only the first order exchange contribution, the 
matching with Eq. (A. 5) allows one to resum the ladder series obtaining 



-Vi 2Vi ^ 
2Vi -V,^ 



(A.9) 



32 



Finally, from the matching of Eq. (A. 5) with all the four diagrams in Figs. (A.l) 
and (A. 2), one gets the full RPA result (direct + exchange): 



^RPA 



0, KY^ 



3Vi 



V 



RPA 



3Vi 
3Vi 



(A.IO) 



B The HF polarization propagator in a broken vacuum 

In this appendix we display how the calculations leading to 11^^ (which enter 
into Eq. (51) for H?^!^'^) can be performed analytically in a quite straightfor- 
ward way following a procedure first introduced in Ref. [25]. Starting from 
the expression for given in Eq. (42a), one introduces the particle-hole 
propagator in the HF approximation, defined as: 

c^F^j: -> . _ o{\k + q\- k-^)eik^ - k) e{k^ - \k + q\)e{k - kt) 

The above expression can be manipulated in the following way: 



k + q]- kp)d{kt - k) e{kp -\k + q\)d{k - k+) 



+ 



9{kp- \k + q\)e{k% - k) 9{kp - \k + q\)e{kt^ - k) 



ejkt - k) - ejkp - \k + q\) 



(B.2) 



being r]^ — sign(a;)?7. Hence, expressing H^^ in terms of the particle-hole 
propagator 

nHF/ [ ^k ^ f dk 9{k+ - k) - 9{k^ - \k + q\) 

J (27r)3 P J (27r)3 cj - u;^^^ + cut + 17]^ 

(B.3) 

— * — * 

and performing the change of variable k + q ^ k within the integral for the 
second contribution, one gets: 



dk 



9{kt - k) 



9{kF - k) 



The energy denominators in the equation above can be written as 



(B.4) 
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' 2m m 

_i_ _ q-k 

a; + a;T - a;- = a; + \ 

^+1 ^ 2m m 



where 



27r2 



- Ao;, 



(B.5b) 
(B.6) 



Note that at equihbrium Vi is related to A;^ through Eq. (6) and one gets to 
the expression (46) of Section 4. 



It is particularly convenient to introduce the scaling variables: 



kp 



muj q 



mAuj 
k$q ■ 



In terms of the above variables we can write: 



(B.7a) 
(B.7b) 



- . ^ kpQ I / A , Q-k 



m 

kpq 
m 



hjp 

Hence 11^^^ can be expressed analytically as follows: 



n^l(g,a;) 



m 1 



(B.8a) 
(B.8b) 



{kpfQ^'Ki^- - M-) - {k^fQ^'\i^+ - AV^+ + q+] 



q (27r)2 

being q± = q/k$. Analogously, for n^|^^ one gets: 



(B.9) 



n;!(g,a;) 



m 1 
7(2^ 



2n{o), 



(B.IO) 

The expressions reported above allow us to write Rell™ (Ren?|^^) and Imll^^ 
(Imll^^) in terms of Legendre polynomials and functions of second kind, Pn 
and Qn, respectively. Indeed it is easy to show that 



ReQ(°)(^) = i[goW-g2W] 



ImQ(°)(V') = -sign(a;)^(l - V'')|[^o(V') - P2 



(B.lla) 
(B.llb) 
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